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1. Introduction 



X 



Despite the very rich physics it offers, the soft gamma-ray sky 
is not fully exploited at present because the telescopes in this 
domain are blinded by intense and complex instrumental back- 
ground induced in their detectors by space environment (cos- 
mic rays, radiations belts, Earth albedo and solar ares, see e.g. 
[Weidenspointner et al., 2005]). Current instruments operating 
in this part of the electromagnetic spectrum do not use focus- 
ing optics. They reconstruct the incidence direction of detected 
events thanks to either an aperture modulation (coded mask) or 
by tracking the multiple (Compton) interactions of photons in 
a sensitive volume. The common point of these two techniques 
is that the signal is collected onto an area which is itself the 
sensitive area. 

To keep improving our knowledge of the violent celestial 
processes responsible of the emission of high-energy photons 
we need to develop more sensitive telescopes. With the existing 
kind of telescopes more sensitive means larger in order to col- 
lect more signal. But the improvement of sensitivity only scales 
with the square root of the collection surface since the instru- 
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Crystals are the elementary constituents of Laue lenses, an emerging technology 
which could allow the realization of a space borne telescope 10 to 100 times more 
sensitive than existing ones in the 100 keV - 1.5 MeV energy range. This study 
addresses the current endeavor to the development of efficient crystals for the 
realization of a Laue lens. 

In the theoretical part 35 candidate-crystals both pure and two-components are 
considered. Their peak reflectivity at 100 keV, 500 keV and 1 MeV is calculated 
assuming they are mosaic crystals. It results that a careful selection of crystals can 
allow a reflectivity above 30% over the whole energy range, and even reaching 
40% in its lower part. 

Experimentally, we concentrated on three different materials (Sii_ A Ge x with gra- 
dient of composition, mosaic Cu and Au) that have been measured both at ESRF 
and ILL using highly-monochromatic beams ranging from 300 keV up to 816 
keV. The aim was to check their homogeneity, quality and angular spread (mosaic - 
ity). These crystals have shown outstanding performance such as reflectivity up 
to 31% at - 600 keV (Au) or 60% at 300 keV (SiGe) and angular spread as low 
as 15 arcsec for Cu, fulfilling very well the requirements for a Laue lens applica- 
tion. Unexpectedly, we also noticed important discrepancies with Darwin's model 
when a crystal is measured using various energies. 

mental background scales with the volume of detectors. This is 
why it appears impossible to make the required sensitivity leap 
of a factor 10-100 with the existing principles of soft gamma- 
ray telescopes. 

A new approach involving the concentration of gamma rays 
has been studied for more than a decade and proved to be fea- 
sible in the ^100 keV - 1.5 MeV domain [von Ballmoos et al., 
2004]. This paper addresses the endeavour that has been ongo- 
ing for the last four years to develop efficient elementary con- 
stituents for such lens: crystals. 



The principle of Laue lenses is described in the next section. 
The third section gives an overview of theories used to model 
diffraction in both mosaic crystals and crystals having curved 
diffracting planes. The fourth section deals with the theoretical 
study that has been undertaken to identify among pure-element 
and two-component crystals the best suited for a use in a Laue 
lens. The fifth section is devoted to the presentation of exper- 
imental characterizations of the three different types of crystal 
of which diffraction properties have been investigated: gradient 
silicon-germanium (SiGe), copper (Cu) and gold (Au). Analy- 
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sis of experimental results and comparison with theoretical pre- 
dictions are given in subsections associated with each crystal 
type. The last subsection emphasizes the discrepancies that we 
have noticed with mosaic crystals between Darwin's model and 
experimental results. The conclusions and perspectives of this 
work are given in the last section. 

2. Principle of a Laue lens 

A Laue lens concentrates gamma-rays using Bragg diffraction 
in the volume of a large number of crystals arranged in concen- 
tric rings and accurately orientated in order to diffract radiation 
coming from infinity towards a common focal point (e.g. [Lund, 
1992]). In the simplest design each ring is composed of identi- 
cal crystals, their axis of symmetry defining the line of sight of 
the lens (see Figure 1). 

Bragg's law links the angle 9 B between the rays direction of 
incidence and the diffraction planes to the diffracted energy E 
through the diffracting planes d-spacing dm- 



2d m sin B = 
<=> 2d HKL sin 9 B = n 



he 
~E 
he 
~E 



(1) 



with h, k, I the Miller index defining the set of diffracting planes 
at work (another notation uses H, K, L prime numbers and n the 
order of diffraction), h the Planck constant and c the velocity 
of light in vacuum. Considering a focal distance /, the mean 
energy diffracted by a ring only depends on its radius r and the 
d-spacing of its constituent crystals [Halloin & Bastie, 2005]: 
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(2) 



We can distinguish two types of lenses. The first one uses 
either various crystalline materials or various reflections (keep- 
ing in mind that higher orders are less efficient) to maintain 
the product dm r constant in every ring so that they all diffract 
the same energy {E\ — E2 in Figure 1). This principle allows 
the realization of a narrow bandpass lens like the prototype 
CLAIRE which used eight Ge reflections in eight rings to focus 
around 170 keV [Halloin et al., 2003; von Ballmos et al., 2004]. 

On the other hand if the same reflection of a given material is 
used on many consecutive rings, the diffracted energy is slightly 
shifted from one ring to the next increasing towards smaller 
radii (E\ > E-i m Figure 1). If the crystals diffract a bandpass 
large enough, the contributions diffracted by neighboring rings 
overlap resulting in a broadband continuous energy coverage. 

Recent projects of Laue lens telescope combine both effects. 
They use several ranges of rings (each range being composed of 
rings using a single reflection) that superimpose to cover effi- 
ciently broad energy bands. One such project, the Gamma Ray 
Imager (GRI) that has been proposed 1 to the European Space 



Agency (ESA), uses mainly Cu 1 1 1, Cu 200, Cu 220 and SiGe 
1 1 1 to create two broad bandpasses ranging over 220 - 650 keV 
and 790 - 910 keV achieving several hundreds of cm 2 of effec- 
tive area [Knodlseder et al., 2007; Barriere et al., 2007]. 

Single perfect crystals only diffract an energy band of a few 
eV wide, given by the Darwin width (see e.g. [Authier, 2001]). 
A Laue lens requires that crystals diffract a larger band. Two 
types of crystals can achieve this: mosaic crystals and crystals 
having curved diffracting planes (CDP crystals). For a GRI-like 
lens crystals optimal angular spread (hereafter called mosaic- 
ity 2 ) has been proven to be around 30 arcsec [Barriere, 2008], a 
figure which should be multiplied by a factor of ^2 - 3 for Laue 
lenses having a shorter focal length (GRI has a focal length of 
100m). 

3. Theory of diffraction in mosaic and CDP crystals 

3.1. Definitions 

We shall introduce two quantities that will be used in the fol- 
lowing. The reflectivity is defined as the ratio of the diffracted 
beam intensity over the incident beam intensity. And the diffrac- 
tion efficiency is defined as the ratio of the diffracted beam 
intensity over the intensity of the transmitted beam when no 
diffraction occurs. 

3.2. Mosaic crystals 

Mosaic crystals are described using Darwin's model as an 
assembly of tiny identical small perfect crystals, the crystal- 
lites, each slightly misaligned with respect to the others accord- 
ing to an angular distribution usually taken as Gaussian [Dar- 
win, 1914, 1922]. Zachariasen [1945] gives the equation of the 
intensity diffracted as a function of the angle of incidence: 



h,mos = h ^ (l - e 



-2aT„ 



(3) 



where Iq is the incident intensity, To is the crystal thickness, \i 
is the linear absorption coefficient, and a can be interpreted as 
the coherent diffusion coefficient that is written as: 



a = W(A6)Q 



(4) 



where A8 is the difference between the actual angle of incidence 
of the beam onto the diffracting planes and Bragg angle. W is 
the distribution function of the crystallites orientation and Q is 
the integrated intensity diffracted by a single perfect crystal per 
unit of thickness. Q is given by the dynamical theory of diffrac- 
tion: 



Qdyn — 



7T 2 dm 



f(A) 



Aq cos 8 b 

where in the Laue symmetric case f(A) is given by: 

B (2A) + |cos20 B |fi o (2A|cos20 J j|) 



M) = 



2A(1 +cos 2 6» B ) 



Bq(2A) 
2A 



(5) 



(6) 



(7) 



1 in June 2007, as an answer to the first announcement of opportunity of the long term plan Cosmic Vision 2015-2025 

2 for the sake of simplicity, we propose to call mosaicity (noted Q.) the full width at half-maximum of angular distribution of diffracting planes, even in the case of 
crystals with curved planes 
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The above approximation (7) can be done for 8 B small, which is 
always valid for energies above 100 keV. Bq is the Bessel func- 
tion of zero order integrated between and 2A and A is defined 

as 



it to 



A cos 6 B ' 



(8) 



in which to is the crystallites thickness. Ao is called the extinc- 
tion length and is defined for the Laue symmetric case (see e.g. 
[Authier, 2001]) as: 



Ao 



7T V c COS 9 B 

r e \ \C\ \F, 



(9) 



hkl 



with Fhu being the structure factor (taking into account the elec- 
trons repartition in space, the crystal lattice and the effect of 
temperature via the so-called Debye-Waller's factor which is 
included in the structure factor), V c the volume of the crystal 
cell, and C the polarization factor. 

We note that the dynamical theory tends towards the kine- 
matical theory when f <C A (which implies f(A) — > 1 in (5)). 
In this case crystals are referred to as ideally imperfect crystals 
and Q is given by Q kin = A f c ^ . 

The distribution function of the crystallites orientation W can 
be expanded as: 



(10) 



where Q. is the full width at half maximum (FWHM) of the 
angular distribution of crystallites (called mosaicity or mosaic 
spread). 

The reflectivity is derived from equation (3), with its peak 
value obtained for A9 = 0: 



^peak 



I h^,s( T o) _ 1 
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(\-e- 2w ^)e^ (11) 



The thickness maximizing the peak reflectivity is derived from 
(11): 



-jnpeak 
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(12) 



Halloin & Bastie [2005] propose a comprehensive version of 
the calculation of the diffracted intensity in a mosaic crys- 
tal. Specifically the equations to compute the structure factor 
(including atomic form factor and Debye- Waller factor) and a 
numerical method to calculate the Bo function. 

3.3. Crystals having curved diffracting planes 

CDP crystals can be obtained in three ways [Smither et al., 
2005]. One is by applying a thermal gradient perpendicular to 
the considered planes of a perfect single crystal. The second 
is by bending elastically a perfect single crystal. This bending 
can be obtained thanks to an external device applying a force 
on the crystal (as it is commonly made for monochromators in 
synchrotron radiation facilities), but also through the deposition 



of a coating on a wafer or by grinding or grooving one face of 
a wafer. The third way is by growing a two-component crystal 
whose composition varies along the crystal growth axis. 

Thermal gradient gives excellent results because it pro- 
duces a very pure spherical curvature, but it requires a signif- 
icant amount of power which is not available onboard a space 
borne observatory. Elastic bending by surface treatment seems 
promising for a Laue lens but it has not yet been investigated. 
Our CDP crystals were obtained by growing composition gradi- 
ent crystals, which produces intrinsically curved planes without 
any mechanical stress. 

Equations describing the diffraction in such crystals are given 
by Malgrange [2002]. They are an extension of the PPK the- 
ory of diffraction in distorted crystals [Penning & Polder, 1961; 
Kato, 1963] for the case of a large and homogeneous curvature. 
In this theory the distortion of diffracting planes is described by 
the strain gradient j3: 



= 



Ao d 2 h.u 



cos 2 9 B dsodsh 



(13) 



which can be written in a simpler way in the case of a uniform 
curvature: 



(14) 



Q. being the FWHM of the angular distribution of planes, or 
mosaicity 2 , To the thickness of the crystal and 6 W half the Dar- 
win width (Darwin width is defined as 2S W ). 

When the strain gradient becomes larger than a critical value 
(3 C — Balibar et al. [1983] showed that a new wavefield 
is created which decreases the intensity in the diffracted wave- 
field. In the case of a uniform curvature of planes, when the con- 
dition j3 > p c is fulfilled the intensity diffracted at the plateau is 
given by Malgrange's formula [Keitel et al., 1999; Malgrange, 
2002]: 
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The a parameter appearing in (16) can be expanded as 



o = — 



7T \p\ _ de/dTAo a a 2 , 
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(15) 
(16) 

(17) 



It is interpreted as the angular variation of the diffracting planes 
orientation over the extinction length, in unit of Darwin width. 

The thickness maximizing the reflectivity can easily be 
obtained under the hypothesis that the curvature of planes c p 
is uniform. In this case c p can be related to the mosaicity and 
the thickness as follows: 



Inserting (18) in (16), one gets 



h,curved = k l-e o <? 



(18) 



(19) 
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Let us define 



M = 



n d hk i 



and N = 



/i£2 

COS On 



(20) 



It follows that the curvature of planes maximizing the reflectiv- 
ity is obtained by solving the equation: 



34 
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(21) 
(22) 



Using (18) we obtain the thickness maximizing the reflectivity 
as a function of energy and mosaicity (for a given crystal and 
reflection): 



To 



nin(i + f) 



M 



(23) 



3.4. Comparison 

There are two main differences between diffraction proper- 
ties of mosaic and CDP crystals. Firstly, the diffraction effi- 
ciency is limited to 50% in the case of mosaic crystals while 
it can reach 100% with CDP crystals. Secondly the diffraction 
profile of a mosaic crystal is close to the Gaussian while it is 
rectangular for a CDP crystal. The spatial distribution of the 
beam diffracted by a crystal being the convolution of its diffrac- 
tion profile by its spatial extent (assuming a uniform beam 
larger than the crystal), the footprint of a CDP crystal onto the 
focal plane does not have the large tails of a Gaussian. In other 
words CDP crystals better concentrate signal than mosaic crys- 
tals. These two points together make a sizable difference in the 
resulting sensitivity of a telescope. In fact, a telescope having 
its lens composed of CDP crystals would be 75 % more sen- 
sitive than the same telescope which lens would be composed 
of mosaic crystals (considering identical materials) [Barriere, 
2008]. 

Hence CDP crystals are substantially more suitable than 
mosaic crystals for the realization of a Laue lens, but in real- 
ity they are much more difficult to produce. Sii_ x Ge x are the 
only ones that we have been able to procure up until now. A 
feasibility study concerning Gei-jSn* has been carried out at 
the Institute of Crystal Growth (IKZ, Berlin, Germany) which 
concluded that it was not feasible. Indeed the maximal solubil- 
ity of Sn in Ge is about 1 % (instead of 100% for Ge in Si) which 
strongly limits the range of the composition gradient. Moreover, 
the growth of such crystal is extremely difficult because of con- 
stitutional under-cooling on the liquid-solid interface, so that 
the 1 % maximum concentration cannot be reached experimen- 
tally. Another option that has not been explored yet is Nii_ x Sn x 
(mentioned in [Smither et al., 1982]). Seen the reflectivity of Ni 
(Figure 3 in the next section) the potential of such a composition 
gradient crystal appears very appealing. 

4. Search for suitable crystals: theoretical study 

Our field of investigation has been limited to pure materials and 
two-component crystals. Crystals composed of more than two 



elements have not been considered because their large crystal 
cell volume strongly decreases their diffracted intensities. We 
will start with the selection of potentially interesting pure mate- 
rials. 

Suitable materials to realize a Laue lens must first and fore- 
most exist at crystalline state at ambient temperature and pres- 
sure conditions without being too reactive in air (spontaneous 
combustion, deep oxidation) nor radioactive and overly toxic. 
Secondly they must diffract efficiently X-ray radiations, which 
is linked to a high electron density and depends also on the 
crystal lattice: the most efficient crystals have either a diamond 
(diam), a face-centered cubic (fee) or a body-centered cubic 
(bec) lattice. Figure 2 shows the electron density of pure ele- 
ments in crystalline state. Half-filled squares represent elements 
fulfilling the above conditions but having a melting point above 
2000 °C and/or being rare and expensive. Filled circles repre- 
sent elements fulfilling the conditions and being more readily 
available than the ones represented as half-filled squares. 

Eighteen out of about one hundred elements are potentially 
interesting, namely: Al, Si, V, Cr, Ni, Cu, Ge, Mo, Rh, Pd, Ag, 
Ba, Ta, W, Ir, Pt, Au and Pb. Among these eighteen elements 
numerous are soft and/or ductile as for instance Cu, Ag, Au, 
Pb. Crystals must be mechanically robust enough to undergo an 
accurate orientation and must not deteriorate during the intense 
vibrations induced by the rocket launch. Soft materials are nev- 
ertheless kept in the list because a slight doping is often enough 
to modify mechanical properties without changing the X-ray 
diffraction properties. 

Another important selection criterion is the availability in 
large quantities with a good homogeneity. Crystals whose 
growth is well mastered are of prime interest. Unfortunately 
applications of pure-element crystal are rather limited, which 
explains why there are almost no industrial patterns to get a 
large quantity of constant-quality ingots, at the exception of Ge 
and Si. On the other hand two-component crystals are used in 
many applications like for instance GaAs, InAs, InP in electron- 
ics, CaF2 for UV lithography, or CdTe for hard X-ray detec- 
tors, etc... Therefore we have added to the eighteen pure ele- 
ments an arbitrary selection (still open to new entries) of two- 
component crystals easily available in industry that may repre- 
sent an interesting perspective for the realization of a Laue lens. 
An extra advantage of two-component crystals is that some of 
them cleave. To have the external faces representative of crys- 
talline planes can be a strong asset for the mounting and orien- 
tation of a large numbers of crystals on a lens structure. 

In order to quantify and compare the diffraction capability of 
these crystals we have calculated their peak reflectivity for three 
different energies covering our domain of interest (100 keV, 500 
keV and 1 MeV) assuming they are mosaic crystals. Darwin's 
model using the dynamical theory (11) has been employed con- 
sidering that the mosaicity is 30 arcsec, the crystallites thick- 
ness is 5 fim and the thickness is calculated to maximize the 
peak reflectivity according to equation (12), within the limits 1 
mm ^ To ^ 25 mm. In each case the most intense reflection is 
considered as stated in Table 1 . 

Figure 3 shows the results. As expected at high energy the 
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crystals having the highest reflectivity are the ones with a high 
mean atomic number (Z) (see Table 1), and conversely at low 
energy the highest reflectivity is produced by crystals having a 
low mean Z. For the low energy side, this effect is due to the 
fact that the minimum thickness allowed in the calculation is 1 
mm, which makes high mean Z crystals very absorbent at 100 
keV. Various crystals both pure or two-component can yield a 
reflectivity above 40% at 100 keV. In the scope of this study the 
best at 100 keV is MgO, which have two advantages: it cleaves 
along (100) planes, and it is very hard (hardness of 6 on Mohs 
scale). However the possibility to find it with a suitable mosaic- 
ity is still an open question. 

Crystals with hexagonal compact (he) lattice have not been 
regarded in this study because this lattice is not as efficient as 
the cubic one. As an example 44RU (he) and 4 5 Rh (fee) have 
almost an identical density of electrons (Figure 2) but Ru is far 
less efficient than Rh at 500 keV and 1 MeV. 

At 500 keV we see that the two-component crystals are 
not any more a good option. Oppositely Cu, Ni, Ag, Rh, and 
Pb exhibit high reflectivities. They are relatively affordable 
and produced with a mosaic structure (private communication), 
which makes them good potential candidates. At 1 MeV, in 
the conditions taken for this calculation, only Ta, W, Ir, Pt, 
Au and Pb are yielding a reflectivity above 25%, Pb being 
the most efficient bearing more than 30%. The group Ta, W 
and Ir have very elevated melting points (2996 K, 3410 K and 
2410 K respectively) causing them to be expensive and difficult 
to get with a constant quality. Pt is very expensive and hence 
not produced in industrial quantities (several kilograms). At the 
moment remains only Pb and Au as viable candidates. Au is 
quite expensive but still affordable and seems to be available 
with a suitable mosaicity, as shown in the experimental results. 

This study highlights that a careful choice of crystals can per- 
mit to cover from 100 keV to up to 1 MeV with a reflectivity 
above 30 %, which is outstanding for this energy domain. It is 
necessary to identify many candidates for the same energy sub- 
range in order to then use only those whose actual measured 
quality is satisfying our requirements. 

5. Experimental results and discussion 

5.1. Facilities and data treatment 

Since 2005 numerous samples have been measured using two 
facilities in Grenoble (France): the GAMS 4 instrument at ILL 
and the beamline ID 15 A at the European Synchrotron Radiation 
Facility (ESRF). 

At ESRF we used two Ge 711 monochromators bent at the 
Rowland circle to get a sharp monochromaticity between about 
280 keV and 600 keV with a fixed exit. Our samples were 
held by a sucking plate specially designed to allow fast sample 
changes without damage. The measurements consist in rocking 
curves (RCs) in Laue geometry both in diffraction and trans- 
mission. RCs are performed one after the other with a single 
high-purity Ge detector which is moved from one beam to the 
other. The beam intensity is monitored by the current of elec- 
trons in the storage ring of the synchrotron. 



The GAMS 4 instrument uses the gamma-ray flux produced 
by neutron capture in a target inserted close to the nuclear reac- 
tor of ILL. We used an erbium target to produce lines among 
which was the one at 815.986 keV (AE/E ~ 10" 6 ). To select 
the line we used a low mosaicity quartz crystal that gave a beam 
divergence of about 2 arcsec. The same type of detector and 
method as in ESRF was used to record the RCs. 

RCs in transmission and diffraction recorded on the same 
area of the crystal samples are put together and normalized 
by the intensity of the transmitted beam when no diffraction 
occurs. Thus complementary peak and valley show directly the 
diffraction efficiency. Reflectivity is then obtained by applying 
the transmission coefficient to take into account the absorption 
through the crystal. 

In the case of mosaic crystals RCs are fitted using Darwin's 
model (as defined in section 3.2) which allows the extraction 
of the mosaicity and of the crystallites size. To be coherent the 
diffraction efficiency values presented in this paper also come 
from the fit of Darwin's model. In the case of CDP crystals, an 
average of the diffraction efficiency is made over the width of 
the plateau. 

5.2. Copper mosaic crystals 

Copper crystals are produced by the monochromator group 
of ILL, where they manage the growth of 8 kg ingots of very 
low mosaicity [Courtois et al. 2005]. The group has extensive 
experience in growing high quality Cu crystals, but our require- 
ments in mosaicity are far below their usual needs. Mosaic- 
ity and homogeneity were the main challenge with Cu mosaic 
crystals. Lately, some pieces of Cu mosaic crystal featuring a 
mosaicity below 1 arcmin have been measured. Three exam- 
ples are shown hereafter (all samples have been measured as-cut 
without surface treatment). 

The first one is the sample 834. 84 (55 x 20 mm 2 cross sec- 
tion and 3 mm thick). It was measured at ESRF in 25 points 
using the 111 reflection in a 299 keV beam. Averaged values 
are reported in the Table 2 and examples of RCs are shown in 
Figure 4 a). As we can see Darwin's model (continuous line) 
fits quite well the data. Note however that the RCs show small 
enlargements at the base most likely due to defect induced by 
the cut. This sample shows a good homogeneity with an average 
mosaicity of 25 arcsec and a standard deviation of 6 arcsec only 
(over 25 spots). Its average diffraction efficiency is 46 %, which 
corresponds to an excellent reflectivity of 34 %. This sample is 
an excellent example of what is needed to build a Laue lens. 

The second sample, 834.(53 (55 x 20 mm 2 cross section and 
9 mm thick), has been measured at ESRF in 25 spots using 
the 111 reflection in a 589 keV beam. The average parame- 
ters extracted from the data are summarized in the Table 2 and 
examples of the RCs are shown in Figure 4 b). Despite instru- 
mental problems that produced an oscillation of the counting 
in transmission geometry we immediately notice that Darwin's 
model fits very well and that these RCs are very narrow for a 
Cu crystal. Indeed the average mosaicity of this sample is 14 
± 6 arcsec which is among the smallest values ever measured 
in a Cu crystal. Initially the thickness was optimized for 850 
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keV which explains that despite a good diffraction efficiency 
the reflectivity is 'only' 26%. 

The last Cu sample, Cu 834.31 (19.2 x 21 mm 2 , 10.3 mm 
thick), has been measured in the 220 reflection at 816 keV 
(GAMS 4) over ten large spots of 2 mm x 1 1 mm covering more 
than half its face. RCs shown in Figure 4 c) are averaged over all 
spots. As indicated in Table 2 the mean mosaicity is 39 arcsec 
and its standard deviation over the ten spots is 6 arcsec which 
indicates a satisfactory homogeneity. 

Disregarding the influence of the crystallites thickness (i.e. 
in the limit of the kinematical theory) the diffraction effi- 
ciency decreases with W (0)Qkin which is proportional to 
d m \F m \ 2 1 E 2 £1. Knowing that d 2 20 < dm and |F 2 2o| 2 < 
|Fin| 2 one understands why this crystal yields a lower diffrac- 
tion efficiency than the two former, equaling 12%. Figure 4 
d) shows the measured peak diffraction efficiency versus the 
mosaicity for the three Cu samples and for comparison theoreti- 
cal values according to the kinematical theory (calculated using 
eq. (11) with a crystallites thickness of 0.1 but disregard- 
ing the absorption term). As can be seen, at 816 keV even an 
ideally imperfect crystal is not efficient with the 220 reflection: 
with a mosaicity of 40 arcsec, we could have expected at most a 
diffraction efficiency of 16.3%(which would have resulted in a 
reflectivity of 9. 1 %). For comparison, in the same conditions the 
111 reflection would lead to a diffraction efficiency of 35.7%, 
making a reflectivity of 19.7%. 

It results from this study that Cu crystals as produced at ILL 
are usable for the realization of a Laue lens. Between energies 
of <~ 300 keV and <~ 600 keV these crystals give optimal results. 
However during our experimental runs we have noticed that cut- 
ting is critical to the final properties of the crystal sample (per- 
formed with an EDM machine at ILL). It can induce defects 
over depths of the order of 1 mm or more which cannot be 
removed by acid etching. A systematic study would have to be 
undertaken to determine the best compromise between cutting 
speed and degradation of the crystals. 

5.3. Sii ^Ge^ concentration gradient crystals 

Sii-xGe^ alloy with x increasing along the crystal growth axis 
(gradient crystal) is produced at IKZ [Abrosimov, 2005]. The 
increase of Ge concentration deforms the silicon lattice result- 
ing in a spherical curvature of the diffracting planes perpendic- 
ular to the growth axis. In our case, crystals have been grown 
along the [111] direction, which results in a spherical curvature 
of the (111) planes. Since the Ge concentration remains low in 
the crystal the curvature is approximatively proportional to the 
gradient of Ge concentration [Smither et al., 2005]: 

VC Ce w e/Rc (24) 

where VC Gc is the gradient of Ge concentration expressed in 
atomic percentage per centimeter, R c the radius of curvature of 
planes (in meters), and e a constant. The assumption that the 
curvature is regular (spherical) permits to relate the gradient of 
Ge concentration to the mosaicity £2 

VC Ge tte£l/T (25) 



with Tq the thickness of the crystal, e equals 25 in Smither's 
paper. 

The first part of the investigation consisted in the character- 
ization of the diffraction properties of several spare samples 
given by IKZ of which SiGe 10.3 is of particular interest. This 
sample (28 mm x 16 mm, 20 mm thick) has been measured 
twice at ESRF using the 111 reflection. The first time we mea- 
sured 15 points on the crystal axis (along the gradient direction) 
using a 297 keV beam. We noticed that the mosaicity increased 
from the low Ge concentration side to the high Ge concentra- 
tion side as expected. The change in curvature is a convenient 
tool to study the accordance with the theory. Among the RCs 
recorded the most outstanding is shown in Figure 5 a). 

This sample has also been measured in 5 points along the 
growth axis using a 489 keV beam. Results of both series of 
measurements have been gathered in Figure 5 b) where the 
peak (or plateau) diffraction efficiency is plotted as a function 
of a (17), a variable that increases both with E and £2. The 
black continuous line shows the theoretical predictions accord- 
ing to (16) using as input a mean Ge concentration of 2.5 at.% 
assumed constant over the 16 mm width of the sample (it inter- 
venes in the calculation of the d-spacing, the structure factor and 
the absorption factor). Accordance with theory seem to appear 
for larger a values when energy increases. This point has been 
noticed on other sample as well, but is not yet fully understood. 

In 2007 three ingots of SiGe were produced in IKZ with the 
objective of checking the relation between growth parameters 
(initial Ge concentration, pulling speed, shape of the crystal) 
and gradient of Ge concentration, which is evaluated through 
the RC width thanks to equation (25). In this series the ingot 
referred to as SiGe 368 has allowed the extraction of six pieces 
of 15 x 15 x 23 mm 3 as represented in Figure 5 c). The pieces 
have been used as -cut without further surface treatment as it has 
been previously established that nothing changed after a deep 
acid etching (0.5 mm on each face). Every piece was measured 
in 13 different spots at ESRF using a 299 keV monochromatic 
beam and the 1 1 1 reflection (RC in transmission and diffraction 
as shown in Figure 5 e)). Figure 5 d) shows the RC width as a 
function of the position along the growth axis, the low Y side 
having the lowest Ge concentration. £2 increases from 30 arc- 
sec up to about 95 arcsec over the 75 mm length of the various 
pieces, but is relatively constant in the two asymmetric pieces 
la and lb. 

Figure 5 f) is a plot of the diffraction efficiency measured 
on each spot as a function of the RC width £2. On this plot the 
continuous line shows the theoretical prediction calculated from 
(16), but this time using an evolving Ge concentration: using the 
measured Ge concentration at the extremity of the samples (the 
low concentration side of samples la and lb and the high con- 
centration side of sample 5), and knowing the corresponding 
solidified fraction g (see Table 3) we refined the value of the 
segregation coefficient to k = 0.42 (assuming it to be constant 
on this part of the crystal) and plugged it into the Scheil-Pfann 
formula: 

C s Ge = kC Ge , {\-g) k - ] (26) 
where C s Ce is the Ge concentration in the solid phase, and Cc e , o 
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is the initial Ge concentration in the melt. Cc e ,o = 0.05 in the 
case of ingot SiGe 368. With this data we can approximate the 
Ge concentration quite well along the 75 mm of the 6 samples 
(see right column of Table 3) and so determine the theoreti- 
cal diffraction efficiency. The agreement with theory is quite 
good especially for the large values of £1, where the curvature of 
planes is larger. We could have plotted the diffraction efficiency 
versus a which would have spanned from 9 to 35 in this case. 
This would show, again, that the accordance is better for larger 
a values than for small ones. Unfortunately it has not been pos- 
sible to measure the SiGe 368 samples at a different energy (we 
have already recorded 156 RC!). It would have permit to check 
if the 'accordance point' shifts towards larger a for a higher 
energy. The agreement between theory and experiments means 
that the curvature of diffracting planes is very regular. We also 
notice that there is no visible difference between symmetrical 
(2,3,4 and 5) and asymmetrical (la and lb) cut elements which 
bodes well for our application since it means that the whole vol- 
ume of the ingot can be utilized. 

The ultimate objective of this development process is to grow 
crystals having a constant gradient of Ge concentration which 
value would produce a curvature optimized for a RC width of 
30 arcsec at 300 keV. All indications point towards this goal 
being achievable. New samples that have been produced in early 
Fall 2008 will hopefully confirm this during a beamtime slot at 
ESRF foreseen in Spring 2009. 

5.4. Prospective for high-reflectivity crystals: Au 

In the context of the search for more efficient crystals for 
high-energy diffraction, a gold crystal (a disc of 10 mm in diam- 
eter and 2 mm thick) has been purchased from the German com- 
pany Mateck. This crystal, which has been polished, was mea- 
sured at four different energies (299, 399, 494 and 588 keV) at 
ESRF during Spring 2008. For each energy five spots (forming 
a square of 5 mm of diagonal plus its center) of 1mm x 1mm 
were probed using the 111 reflection. Results have been com- 
pounded in Table 4 and a typical RC of this crystal is shown in 
Figure 4 e). 

Results are very positive; the crystal seems homogenous with 
a mosaicity ranging from 16 arcsec at 588 keV to 30 arcsec at 
300 keV and achieve an excellent diffraction efficiency. This 
high-Z material was the first of a new wave of investigations. 
A less successful attempt was also conducted with a Pt crys- 
tal (from Mateck) which showed multiple grains and an over- 
all mosaicity larger than 1 degree. Recently several Ag and Rh 
crystals have been purchased (from Mateck) and will be char- 
acterized soon. 

5.5. Limitation of Darwin's model 

In the previous example (Table 4) the crystallites thickness 
seems to increase as the inferred mosaicity decreases. This 
apparently non-physical behaviour has been observed in other 
measurements of mosaic crystals [Barriere, 2008] as well and 
is not explained by Darwins model. As shown in Figure 4 f), it 
is expected that the FWHM of the RCs decreases when energy 
increases. It is due to the fact that the extinction length (9) is 



proportional to the energy: at low energy the extinction length is 
smaller than crystallites, which implies that the diffraction peak 
is cut due to primary extinction in each crystallite. The peak 
being cut, it results in a larger FWHM. Hence the fit of data 
with Darwin's model does not require to decrease the mosaicity 
to fit this FWHM falloff. But the observed effect is much more 
important than what can be modeled by dynamical effects in a 
crystal described by Darwin's model. 

Moreover the apparent increase of crystallites size is hardly 
explained in the same way: if the apparent size of crystallites 
was driven by the extinction length the value of A (defined 
in eq. 8) should remain constant with respect to the energy 
of measurement, which is not the case. Our interpretation of 
these phenomena is that the crystallites in a mosaic crystal are 
not all identical but have a size distribution. Hence one empha- 
sizes either the small crystallites when the beam energy is low, 
because large crystallites have a lower contribution due to pri- 
mary extinction. On the other hand one emphasizes the large 
crystallites at high energy when the longer extinction length 
bring the crystallites in a regime where they diffract propor- 
tionally to their size (i.e. a regime described by the kinematical 
theory). In this idea, the observed decrease of mosaicity can 
be explained if one assumes that the smaller crystallites have a 
correspondingly larger angular distribution than the larger crys- 
tallites, which seems a reasonable hypothesis. 

A consequence of the increase of crystallites size when the 
energy of measurement increases is that diffraction keep in the 
regime of dynamical theory instead of entering the regime of 
kinematical theory as we could have expected seen the increase 
of extinction length. It means that crystals are less efficient that 
what we could have expected from their parameters deduced 
from low-energy measurements. Accurate performance esti- 
mates of Laue lenses require a reliable modeling of intensities 
diffracted by crystals. Since the energy dependance of the crys- 
tal parameters is not yet modeled, it is not possible to extrap- 
olate accurately the performance of a crystal from one energy 
to another. There are two possibilities to solve this problem, 
either we try to refine Darwin's model integrating a dispersion 
of crystallites size or we can try to determine the 'evolution' of 
the crystallites size that we have to input in Darwin's model to 
stick with experimental results. In the present dataset the evo- 
lution of A with energy is very well fitted by a second degree 
polynomial curve (R 2 = 0.9994) with the parameters: 

A = 3.759 - 7.650xl0~ 3 E + 6.382xl0~ 6 E 2 . (27) 

But with only one crystal measured on one spot at four ener- 
gies it is hard to derive conclusions. Is this evolution of A the 
same for any crystal? Is this equation still valid if we extend the 
energy range of investigation? Either to refine Darwin's model 
or to model the evolution of A, in both cases more measure- 
ments are required. 

6. Conclusions - Perspectives 

The last iteration of Sii-xGe^ production have given very sat- 
isfactory results, paving the road towards the final phase of 
our development program: the production of constant-gradient 
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ingots allowing the extraction of homogeneous pieces having 
a 30 arcsec bandpass and optimized for 300 keV. Cu crystals 
have also shown that our requirements are attainable with pieces 
having a mosaicity as low as 15 arcsec, and reflectivity in accor- 
dance with theoretical predictions. 

A new phase in these investigations started with the charac- 
terization of a gold sample. This 2 mm thick sample has kept 
its promises with an excellent reflectivity near 600 keV. A Laue 
lens such as the one designed for the GRI mission could bene- 
fit dramatically from the enablement of high-Z crystals such as 
gold. In the coming months Ag and Rh crystal samples will be 
measured as well, hopefully enlarging our portfolio of usable 
materials, which is of prime importance for the design of effi- 
cient lenses. 

More generally the identification of efficient crystals opens 
the way for additional experimental tests. Two-component crys- 
tals in particular are of great interest since their growth has 
already been developed for other applications. High-Z pure ele- 
ments are quite rare and expensive, but are mandatory when 
building an efficient lens covering energies higher than 600 
keV. Alternatively, crystals having intrinsically curved diffrac- 
tion planes can enhance dramatically the overall performance of 
Laue lenses. This field, which is largely unexplored so far, could 
also find some interesting application in high-energy monochro- 
mators (X-rays and neutrons) a fact that further strengthens the 
drive for initiating a development process. 
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Table 1 

Miscellaneous information about the crystals considered in this study. 



Crystal name 


Cu 834.54 


Cu 834.53 


Cu 834.31 


Reflection 


111 


111 


220 


Thickness (mm) 


3.0 


9.0 


10.3 


Energy (keV) 


299 


589 


816 


Mosaicity (arcsec) 


25 ±6 


14 ±6 


39 ±6 


Peak diffr. efficiency 


0.46 ± 0.02 


0.47 ± 0.04 


0.12 ±0.02 


Peak reflectivity 


0.34 ± 0.02 


0.26 ± 0.02 


0.06 ± 0.01 


Size of crystallites (p,m) 


60 ± 15 


129 ± 43 


210 ±31 


A — 7T/o/Ao 


1.27 ± 0.32 


1.38 ±0.46 


1.17 ±0.17 



Table 2 

Summary of results of the Cu crystals investigated. Mosaicity, diffraction effi- 
ciency and crystallites thickness are obtained from Darwin's model fits. Values 
are averaged over all measurement spots, and their uncertainty is the standard 
deviation. The bottom row allow the crystallites thickness to be related to the 
extinction length. 



Y(mm) 


g 


Coe,exp (at.%) 


C Ge ,calc (at.%) 





0.192 


2.40 ±0.15 


2.38 


15 


0.347 




2.69 


30 


0.484 




3.08 


45 


0.596 




3.55 


60 


0.692 




4.16 


75 


0.777 


5.05 ± 0.20 


5.01 



Table 3 

Position along the growth axis (Y), solidified fraction (g), measured Ge con- 
centration (Cce. ex P ) and calculated Ge concentration (Co e , calc) in crystal SiGe 
368. 



1 [Cotts & Anderson, 1981] 

2 [Begg, 1976] 

3 [Palchoudhuri & Bichile, 1989] 

4 [Kittel, 1970] 

5 [Kumara Swamy et al., 1996] 

6 [Freer, 1981] 

7 [Matsuo Kagaya & Soma, 1986] 

8 [Stevenson, 1994] 

9 [Fawcett, 1969] 
10 [Bashir, 1988] 
"[Racek, 1973] 
12 [ZubIk, 1976] 
13 [Mavroides & Kolesar, 1964] 
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Energy (keV) 


299 


399 


494 


588 


Mosaicity (arcsec) 


30 ±5 


24 ±3 


24 ±2 


18 ±3 


Peak diffraction efficiency 


0.47 ± 0.02 


0.47 ± 0.01 


0.46 ± 0.02 


0.45 ± 0.02 


Peak reflectivity 


0.12 ±0.01 


0.22 ± 0.01 


0.26 ± 0.01 


0.29 ± 0.02 


Size of crystallites (p,m) 


46 ± 3 


52 ±3 


57 ±4 


65 ±6 


A — n fo/Ao 


2.04 ±0.13 


1.73 ± 0.09 


1.53 ± 0.11 


1.47 ± 0.13 



Table 4 

Summary of the results obtained with gold crystal Au_l using the 111 reflection. Mosaicity, diffraction efficiency and crystallites size are obtained from Darwin's 
model fits. Values are averaged over five measurement spots (1mm x 1mm), and their uncertainty is the standard deviation. Bottom row allows comparison of the 
crystallites size with the extinction length. 
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Figure 3 

Peak reflectivity of various potentially intresting crystals, sorted by increasing 
mean atomic number (Z) (and density in case of equal mean Z). The calcul 
assumes mosaic crystals with a mosaicity of 30 arcsec, a mean crystallites size 
of 5 fj,m and a thickness optimized to maximize the peak reflectivity within the 
limits 1 mm ^ Tq ^ 25 mm. For each crystal, the reflection considered is given 
in Table 1 . 
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Figure 2 

Electron density in crystalline matter. 
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Figure 4 

a) RCs recorded on the crystal Cu 834.54.3 using a beam of 299 keV and the 111 reflection. The crystal is 3 mm thick, surface has been left as cut which may be 
the reason of the tails around the peak. The mosaicity in this point is 29 arcsec and the diffraction efficiency reaches 47% yielding a reflectivity of 35% taking into 
account the transmission coefficient of 0.75. b) RCs recorded on the crystal Cu834.<53.2 using a beam of 589 keV and the 111 reflection. Crystal thickness equals 9 
mm, which makes a transmission coefficient of 0.55. Mosaicity in this point is 14 arcsec. The diffraction efficiency reaches 50%, yielding a reflectivity of 27.5%. 
c) Averaged RCs over ten measurement spots on crystal Cu 834.31 using the 220 reflection and a 816 keV beam, d) Peak diffraction efficiencies measured on the 
three Cu crystals 834.31 (averaged over the ten spots), 834.53 and 834.54 as a function of mosaicity. Experimental results are compared to Darwin's model using 
the kinematical theory, e) RCs recorded on 1 1 1 reflection of the gold sample Au_l with a beam energy of 588 keV. The estimated mosaicity is 16 arcsec. Taking into 
account the absorption through the thickness of 2 mm a reflectivity of 31% is found. Continous line indicates a fit to Darwin's model, f) Simulated rocking curves 
using Darwin's model of an Au crystal (111 reflection) of 2 mm thick, 30 arcsec of mosaicity and having crystallites of 40 fim of thickness. 
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Figure 5 

a) One of the most striking RC recorded on the sample SiGe 10.3 from 111 reflection. The thickness of the sample is 20 mm, and the beam energy is 297 keV. 
The mosaicity equals 12 arcsec and the reflectivity 60%. b) Diffraction efficiency measured on crystal SiGe 10.3 at 297 keV (red) and at 489 keV (blue), along the 
growth axis. The continuous line shows the theoretical prediction assuming a constant Ge concentration of 2.5 %. c) Ingot 368 of Sii_ A -Ge x , x varying along the 
growth axis (x increase towards the right in the picture), produced at IKZ. Six pieces each measuring 15 x 15 x 23 mm 3 were extracted (as shown in the figure). The 
dashed line represents the (111) spherical planes, d) Mosaicity as a function of crystal position along the growth axis in the six pieces extracted from ingot SiGe 
368. Measurements have been performed at ESRF using the 111 reflection and a 299 keV beam, e) A very good example of a measured RC from the 111 reflection 
of sample SiGe 368 lb: the rectangular shape of the curves is very close to what is expected in the ideal case, showing the regularity of the curvature of diffracting 
planes. The thickness of the sample is 23 mm and the beam energy was 299 keV. f) Peak diffraction efficiency measured on the six pieces extracted from ingot SiGe 
368 as a function of mosaicity and compared to theoretical predictions (continuous line). 



